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We investigate the spectral dimension obtained from non-local continuum d’Alembertians derived 
from causal sets. We find a universal dimensional reduction to 2 dimensions, in all dimensions. 

We conclude by discussing the validity and relevance of our results within the broader context of 
quantum field theories based on these nonlocal dynamics. 

PACS numbers: 


I. INTRODUCTION 

In recent years, the study of the spectral dimension as 
a tool to investigate the small scale structure of space- 
time has attracted much attention. Indeed, analyses per¬ 
formed in causal dynamical triangulations [T], asymptot¬ 
ically safe quantum gravity [2j, loop quantum gravity [3], 
super-renormalizable gravity [Ml, noncommutative ge¬ 
ometry mm and in cosmology [5], all point towards a 
running of the spectral dimension at short scales [10]. In 
particular, dimensional reduction appears to be ubiqui¬ 
tous in quantum gravity (see mi)- 

Although several approaches to quantum gravity share 
the property of dimensional reduction on small scales 
mi, the causal set approach is an outlier. Indeed, results 
by Eichhorn and Mizera [12] (EM from here on) show 
that the spectral dimension of a causal set increases at 
small scales. Their analysis relies on computing the spec¬ 
tral dimension from the return probability, P , of a ran¬ 
dom walker on a causal set with the spectral dimension 1 
defined as cLem = — 2d In P(s)/d Ins. They then argue 
that the increasing behaviour of cIem at short scales is 
due to the radical nonlocality present in causal sets mi- 
This exact same non-locality which, as we will see shortly, 
arises because of the interplay between Lorentz invari¬ 
ance and discreteness, also appears in the construction of 
d’Alembertians governing the dynamics of scalar fields on 
causal sets [14]. An obvious question then is whether the 
spectral dimension computed from such d’Alembertians 
using the usual heat kernel techniques, corroborates the 
EM result. 

The rest of the paper is organised as follows. In Sec¬ 
tion 2 we briefly introduce causal sets, the notion of 
sprinklings and how the non-locality referred to above 
arises within this context. We then introduce the fam¬ 
ily of d’Alembertians defined in mi in both position and 
Fourier space, and discuss some of their properties. Sec¬ 


1 They also defined a “causal” spectral dimension using the meet¬ 
ing probability of two random walkers, Pm, both moving forward 
in time, which give similar results. 


tion 3 is devoted to computing the spectral dimension of 
Minkowski spacetime. In particular, we determine both 
the short and large scale behaviour of the spectral di¬ 
mension in all dimensions analytically, and for particular 
choices of nonlocal d’Alembertians in 2, 3 and 4 dimen¬ 
sions we numerically compute the full dependence of the 
spectral dimension on the diffusion parameter. In Sec¬ 
tion 4 we discuss our results, and in Section 5 we spell 
out some conclusions. 


II. CAUSAL SETS, NON-LOCALITY AND 
D’ALEMBERTIANS 

Before we define the d’Alembertians that will be used 
in calculating the spectral dimension we briefly review 
how the non-locality referred to in the introduction 
arises. 


A. Causal Set Nonlocality 

Recall first that a causal set (causet) is a locally finite 
partially ordered set, where the partial order is taken 
to underly the macroscopic causal order of spacetime 
events mi- In order to unite fundamental discreteness 
with Lorentz invariance a kinematic randomness must be 
introduced. This is done via the notion of a sprinkling, 
i.e. a Poisson process of selecting points in a spacetime 
M uniformly at random at some density p = 1 /l d , such 
that the expected number of points in a spacetime re¬ 
gion of volume V is pV. This process generates a causet, 
C, the elements of which are the sprinkled points and 
the order relation of which is the spacetime’s causal rela¬ 
tions restricted to the sprinkled points. The randomness 
of the sprinkling is key in ensuring that the causet thus 
generated is Lorentz invariant Em, and it leads to the 
following radical form of non-locality. Consider a causet 
C that is a sprinkling of d dimensional Minkowski space- 
time, and let p £ C. Then the nearest neighbours to, say 
the past of, p will be the set of points q past-linked to p, 
i.e. {q A p | $ r s.t. q A r A p}. The locus of such points 
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will roughly comprise the space-like hyperboloid, E, to 
the past of p defined by E = {q € J~{p) | r(p, g) = ^}. 
Hence there will be an infinite number of such points. 
Precisely this form of non-locality is evident in the defi¬ 
nition of d’Alembertian operator on a causal set, defined 
by constructing a finite difference equation in which lin¬ 
ear combinations of the value of the field at neighbouring 
points are taken. Since the number of nearest neighbours, 
next nearest neighbours, etc., is infinite, the correspond¬ 
ing expression looks highly non-local (see [13] and [14]). 


in ng, to include an infinite family of nonlocal 
d’Alembertians for any dimension d. The nonlocal 
d’Alembertians Dp d \ labelled by a non-locality scale p 
originally taken to coincide with the fundamental dis¬ 
creteness scale of causal sets, but need not be — which 
we take to be a free parameter of the theory, are given 

by 


B. Nonlocal d’Alembertians 

Recently, Aslanbeigi et. al ([12]) generalised the orig¬ 
inal constructions of nonlocal d’Alembertians in m 


= P 1 (a^ d U(x) + p/3^ J2 f d d y {pV{X \ e ~» v ^cf>(y) 


(1) 


where Nd is a dimension dependent positive integer, 
p = 1 /l d and the coefficients a^ d \ fd^ and Cn can be 
found in Equations (12)-(15) of [Tjj3|. It is straightforward 
to recover the original results of [14] in d = 2,4 as the 


simplest cases of the analysis reported in m- 

The momentum space representation of these expres¬ 
sions, g^p\k 2 ), were also computed in [IS] : 


( „ d Nd 

g W(k 2 ) = p 2 ' d a W+pW2(2n) d / 2 ~'Z~^ ^ ° r 


E- 

n—0 


n\ 


Id 


dss d{ - n+1/2) e-' 1dsd K !k _ 1 (Z l / 2 s)^ , 


( 2 ) 


where 


r_k-k (f) 2 /oN 

P 2 ' d ’ ld dT(^±l)’ () 

and K u is the modified Bessel function of the second 
kind, possessing a cut along the negative real axis and 
for which we assume the principal value. 

The functions |2]) simplify considerably in the limits 
where k 2 <C p 2 ^ d and k 2 3> p 2 / d , he. the infrared (IR) 
and ultraviolet (UV) limit respectively. By construc¬ 
tion the IR behaviour is the same for all operators, i.e. 
gp d \k 2 ) —> — k 2 as k 2 —*■ 0. The UV limit on the other 
hand is dimension dependent, and can be shown to be 

g (d \k 2 ) -»■ ap 2 ' d + bp-* +1 {k 2 )~ d ' 2 + ..., (4) 

as k —> oo, where a and b are dimension dependent con¬ 
stants. Note that since g^ goes to a constant in the UV, 


the Green function ([T]) possesses a delta function diver¬ 
gence in the coincidence limit. This divergence can be 
regularised by subtracting the constant ( ap 2 / d )~ 1 from 
the momentum space Green function m- Inverting back 
gives a regularised momentum space d’Alembertian 
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(d) , = 

reg ' 


ap 2/d g d 

ap 2 / d — g 


d ’ 
P 


(5) 


which now has the following UV behaviour 

( 6 ) 

Note that this regularisation procedure is manifestly 
Lorentzian and is physically motivated by the underlying 
theory being a theory on the discrete causal set, where a 
coincidence limit does not exist. 

The operators we started with <0 are retarded Lorentz 
invariant operators, and their Laplace transforms are 
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therefore defined in the limit 5s (k°) —» 0 + [151 1201121j . 
i.e. in the upper half complex fc°-plane. Their inverse 
Fourier transforms therefore yield the (unique) retarded 
Green functions to 0 and are given by the integral 

, f dk° f d^ 1 k e ik < x ~^ 
X ' y) -J TR ^J (27T)-- 1 fl(«0(jfe2)’ 

where is a contour running from —oo to oo in the 
upper half complex k° plane such that all singularities of 
g _1 lie below the contour, see Fig. [l] 

In the next section, in order to compute the spectral di¬ 
mension using heat-kernel techniques, we will need to to 
Wick rotate the d’Alembertian, g regi or equivalently its 
retarded propagator. However such a Wick rotation can¬ 
not be performed on the retarded propagator because the 
contour, T#, would cross singularities. To avoid this issue 
one must use the Feynman propagator whose contour can 
be Wick rotated without crossing any singularities (see 
Fig. [I]). To define this propagator we first analytically 
continue g to the whole complex plane and then define 
the Feynman propagator to be 0 with g reg replaced by 
its analytically continued version and replaced by Tj?. 
Although this propagator is not a Green function of the 
original retarded d’Alembertians, which indeed only ad¬ 
mit unique retarded propagators, two independent stud¬ 
ies of free scalar QFTs based on such nonlocal dynamics 
suggest that it is the correct Feynman propagator for 
such theories (STJ [22]. From here on we drop the explicit 
p dependence. 




FIG. 1: A schematic diagram of the singular structure of 
l/g(k 2 ) in the complex k° plane, including two contours, F_r 
and r p, which define the retarded and Feynman propagators 
respectively. Note that before we are able to define the latter 
we must first analytically continue g{k 2 ) to the lower half the 
plane. The contour F^ will clearly encounter an obstruction 
under a Wick rotation, k° —> —ik°, while Ff won’t. 


III. SPECTRAL DIMENSION 


A. Analytical results 


We first recall that the spectral dimension is defined 
as 


d s 


ain(PQQ) 
5 In (s) 


( 8 ) 


where P(s) := Tr(P(s; x , y)) is the return probability for 
a particle, the diffusion process of which is dictated by 
the Laplacian of interest via the equation 

■^P(s;x,y) + AP(s;x,y) =0. (9) 

Therefore, given a modified d’Alembertian, one has to 
first find the corresponding Laplacian in Euclidean sig¬ 
nature, and then consider the heat equation, on the man¬ 
ifold of interest, determined by the Laplacian at hand 
[251 . Equation 0 can be interpreted as describing dif¬ 
fusion of a fictitious particle on the euclideanised space- 
time, the solutions of which are given by the heat kernel. 
The latter, denoted as P(s;x,y), physically represents 
the probability density of diffusion from x to y in time 
s. The trace of P then gives the return probability for 
the diffusing particle. In flat d-dimensional Minkowski 
spacetime with standard local d’Alembertian for exam¬ 
ple, the spectral dimension coincides with the Hausdorff 
dimension d of the spacetime, for all diffusion times. 

Consider now a free massless scalar field in d- 
dinrensional Minkowski spacetime, with dynamics de¬ 
fined by a member of the regularised nonlocal 
d’Alembertians, as introduced in Section[nJ namely glt\. 
The trace of the heat kernel, P(s), is given by 


P(s) = J J^^ eS9 ™ a = Cd J Q dkk<i 1 e S3 ™ 9 , ( 10 ) 

where Cd is a dimension dependent constant, and we have 
analytically continued to Euclidean signature. 

The IR behaviour of gicg —> — k 2 , as k 2 —> 0, ensures 
that for large diffusion times the spectral dimension will 
flow to the value of the Hausdorff dimension d. In the 
UV instead we have that gi^g « k d , so that the trace of 
the heat kernel can be written as 

/•OO i 

/ dkk d - 1 e~ kds = — (e~ xds ), ( 11 ) 

Jx sd 

up to irrelevant numerical factors. In the above intergral 
we have introduced an IR cutoff A, which has been chosen 
large enough that the UV approximation for gi d \ is still 
valid, i.e. A Z -1 . Substituting this back into 0, we 
find 


In this section we provide both numerical and an¬ 
alytical results for the spectral dimension computed 
from the regularised Laplace transform of the non-local 
d’Alembertian gl%- The numerical analysis will be re¬ 
stricted to the minimal cases in 2, 3 and 4 spacetime 
dimensions d. 


d u s v = 2 + 2X d s -a 2, ass—^0. (12) 

This shows that the improved UV behaviour of the reg¬ 
ularised nonlocal Green functions leads to dimensional 
reduction to d s = 2 for all spacetime dimensions consid¬ 
ered. Furthermore, the linear behaviour in s with large 
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coefficient X d is in accordance with the numerical evi¬ 
dence provided in Figure [2j It is important to note that 
(12) is only valid for small values (with respect to the 
non-locality scale l) of the diffusion parameter s <C l. 

Finally, we should pay attention to the fact that the 
universal dimensional reduction to d s = 2, crucially re¬ 
lies on the regularisation of the Green functions. One 
might therefore ask what happens should one decide not 
to perform the regularisation? The first hint that some¬ 
thing will go wrong with the computation of d s , arises 
from the fact that g^ d \k 2 ) —> constant as k 2 —> 0, so 
that the integral in ( fIo| ) clearly diverges. To analyse the 
small scale behaviour of d s we must therefore introduce 
both an IR cutoff A and a UV cutoff A. Then 


r A 

Pis') ~ lim / dk k d ~ 1 e~ CdS 
A—>-0, A—Kx) J ^ 

A d _ \d 

= lim e~ CdS ---, (13) 

A—>-0, A.—^oo Cl 


where we used the fact that in the UV, g ^ goes like 
a dimensionful and dimension dependent constant, Cd, 
which therefore dominates the integral. Substituting this 
into ([8 ) and ignoring issues about the order of limits, 
we find that d s = 2c^s. This result is in accordance 
with the numerical simulations we have performed, and 


clearly leads to a spectral dimension that does not make 
much physical sense. Indeed, it does not asymptote the 
Hausdorff dimension d for s l, but actually diverges. 

One could have argued from the beginning that with¬ 
out regularisation the spectral dimension would be mean¬ 
ingless, and that it would possess a linear dependence in 
s for all s. For instance, consider the spectral dimension 
in a spacetime with d = 2. Then, 


-2s 


f 0 °° dz g 2 {z)e^ s 
/ 0 °°dz e®W« 


= 4 ps, 


(14) 


where z = k 2 . To obtain (14) we split both the numer¬ 


ator and denominator by introducing an IR cutoff L big 
enough to approximate the integrand with its UV con¬ 
stant behaviour and use the fact that / 0 °° dz e 92 ^ s di¬ 
verges. The factor of 4 obtained in the final expression is 
specific to the minimal non-local d’Alembertian in d = 2. 


B. Numerical results 

We numerically calculate the spectral dimension in 2, 
3 and 4 dimensions using the regularised minimal non¬ 
local d’Alembertians. These are given in terms of the 
unregularised d’Alembertians m- 


2 h (2) 

g^ik 2 )=a^p + 2pY / \ 


n =0 


< __ » n poo 

/7M / dCe n+1 e~^ e K 0 


2 ^3^ 

g {3 \k 2 ) = a^p 2 ' 3 + 2(27r) 1 / 2 /3 5 / 6 (fc 2 ) -1 / 4 £ ^ (A)" f°° d£ , 

n=0 U ' 0 

9 ™<t 2 )=« <4i ^+4*0=) y i 

in d = 2,3 and 4 respectively, where K n are modified Bessel functions of second type and 

a (2) = -2, b™ = 4, b i 2) = -8, = 4, 


(15) 

(16) 
(17) 


a< 3 > = - 


1 1 

( n ) 

V t <3)_ 1 , 

( 7T A 

r (5/3) 1 


1 ’ 0 r(5/3) 1 

UA/ 


2/3 


b 


(3) 

1 


27 (* V /3 1 , ( 4) _ 9 (* A 3 i 

16 \2y/3j r(2/3)’ 2 8\V6j r(2/3)' 


a (4) - A 5(4) - A 
“ A 0 “ A 



36 r(4) _ A 7,(4) _ 32 

A 2 " A’ 3 “ A’ 


as obtained via equation ([ 5 ]). The results given in Fig- approaches 2 in the limit s —> 0, increases to a maxi- 
ure [2] show that in d = 2,3,4 the spectral dimension 
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mum greater than the large scale Hausdorff dimension 
when s ~ 0(1), and then decays back to the Hausdorff 
dimension as s —> oo. 2 


d. 





FIG. 2: From top to bottom we have the spectral dimension 
d s as a function of diffusion time s for p = 1 for d = 2,3 and 
4 respectively. The dashed line represents the value of the 
Hausdorff dimension. In the 2d plot the spectral dimension 
can clearly be seen to interpolate between d a = 2 at short 
scales and d s = 2 at large scales, with a maximum occurring 
for s s» p 1 ' 2 . The large scale asymptotic behaviour and the 
maximum are not as clear in the 3 and 4 dimensional cases 
due to poorer numerics in d > 2. Nonetheless one can see 
the short scale limit to d a = 2, a maximum for s w 10 
and a large scale limit which is very slowly asymptoting the 
Hausdorff dimension. These numerical results corroborate the 
analytical analysis of Section |TTT| and provide evidence for the 
behaviour of d s interpolating between short and large scales 


2 Note that in d = 3,4 the maximum occurs on scales s ~ 10Z 
rather than Z and in the s —> oo the decay back to the Hausdorff 
dimension is not as marked as in 2d. We believe that these minor 
differences are due to the higher degree of numerical approxima¬ 
tion present in their analysis. 


IV. DISCUSSION 

We have studied the spectral dimension arising from 
non-local d’Alembertians derived from causal-set the¬ 
ory. We have shown that after regularising these 
d’Alembertians by removing an unphysical 3 divergence 
in the coincidence limit, dimensional reduction is present. 
In particular the spectral dimension goes to 2 at short 
scales in every dimension. The small scale dimensional 
reduction can be seen to arise from the improved UV 
behaviour of the propagator for such non-local theories. 
We have also provided numerical evidence for d s as a 
function of diffusion time s for the minimal regularised 
d’Alembertians in d = 2, 3 and 4, confirming our analyti¬ 
cal estimates in the limits s <C l and s l. These simula¬ 
tions show that in all three cases considered the spectral 
dimension possesses a maximum on scales of order l, and 
asymptotes the Hausdorff dimension from above. This is 
an indication of super-diffusive behaviour on intermedi¬ 
ate scales — where effects due to the non-locality scale 
start to become relevant — a feature similar to that found 
in |24| in the context of causal dynamical triangulations. 

V. CONCLUSIONS 

To conclude with we would like to speculate on the 
significance and validity of these results. Recall first that 
the UV behaviour —► (fc 2 ) d / 2 was key in ensuring 

that d 8 (s) —I 2 as s —> 0 in all dimensions. Now, it was 
shown in [IS] and [22], where the quantum theory of free 
scalar fields with nonlocal dynamics was studied, that 
these operators lead to Wightman functions of the form 

W(x,y) =J^dp 2 p(p 2 )j d ^ [ e(k 0 )5(k 2 + p 2 )e ik ^~y), 

where p(—k 2 ) = lm(gO)(k 2 ))/\gO)(k 2 )\ 2 , a form reminis¬ 
cent of the Kallen-Lehman representation of the Wight¬ 
man function in interacting QFTs. 4 Unlike the Kallen- 
Lehman representation of local interacting QFTs, how¬ 
ever, the spectral function p arising from the causet de¬ 
rived nonlocal d’Alembertians defined in [TS] is not a pos¬ 
itive function in general. In fact, evidence suggests that 
for d > 2 all operators of this kind will lead to a non¬ 
positive spectral functions p. This implies the existence 
of negative norm states in the quantum description of 
these theories. 

As was explicitly stated by Weinberg (see Ref. [25] Sec¬ 
tion 10.7, p.460), a consequence of the Kallen-Lehman 


3 From the point of view of the fundamental discrete theory, in 
which the coincidence limit is ill-defined. 

4 We have ignored contributions coming from complex mass poles 
in the propagator since they are not relevant for the current 
discussion and their inclusion does not change the conclusion of 
this argument. 
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spectral representation, and the positivity of the spectral 
function p(p 2 ) (which in standard local QFT is given 
by the sum of squares of matrix elements of the field- 
observable), is that the propagator cannot vanish for 
\k 2 \ -A oo faster than the bare propagator 1/k 2 . Since 
the Wightman function of nonlocal theories of the kind 
studied in this paper can be written a la Kallen-Lehman, 
a natural question is whether the improved UV behaviour 
of the propagator, and therefore the small scale dimen¬ 
sional reduction in the spectral dimension, are inextrica¬ 
bly linked to negative norm states in the quantum theory. 
If this were the case, it would raise doubts as to the phys¬ 
ical relevance of our result in the regime s <C l. 

Knowing about the issue of non-positive spec¬ 
tral functions associated to causet derived nonlocal 
d’Alembertians, Aslanbeigi and Saravani studied QFTs 
of generalised versions of the operators in [15], where 
they imposed the positivity of p(g 2 ) from the beginning 
EH- Hence, from Weinberg’s argument one would ex¬ 
pect that their newly defined gl%(k 2 ) goes like k 2 as 
k 2 -A oo, which is indeed the case. Because of this the 
spectral dimension does not reduce to 2 at small scales 
but rather starts off at 4, increases, and then asymptotes 
back to 4 from above at large scales (much like the 2d case 
considered in Section III BI. Although interesting in its 
own right, the specific theory they consider does not fall 
within the family of d’Alembertians used in this paper, 
and as such cannot be (at least trivially) seen as contin¬ 
uum approximations to the more fundamental operators 
living on an underlying causal set, in the way Sorkin had 
originally constructed them. 

A possible way out of these issues arising in the quan¬ 
tum theory is to take the non-locality scale l to corre¬ 
spond to the causet fundamental discreteness scale. In 


this case one would expect that on scales s ~ l the dif¬ 
fusing particle would cease to see a continuum space- 
time, and would start feeling the underlying discreteness 
of the causal set. Therefore, at these scales the analy¬ 
sis performed in this paper would cease to be applicable. 
In order to study the behaviour of the spectral dimen¬ 
sion on scales s < l, one would instead have to resort 
to methods to be deployed directly to the causal set it¬ 
self, e.g. Eichhorn et al.’s method [TJ. We conclude 
by noting that our analysis shows an increasing d s (s) as 
s —> Z + , which then only starts to decrease after it has 
gone through s ~ 0(1). This might be a hint of the exis¬ 
tence of a universal description of the spectral dimension 
which interpolates between our results and the EM spec¬ 
tral dimension. 


Note added 


After this work was completed the article [SB] by 
S. Carlip appeared in the electronic archives with sim¬ 
ilar results. 
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